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Abstract 

The dynamic dipole polarizabilities for the four lowest triplet states {2^S, 3^S, 2^P and 3^P) 
of helium are calculated using the B-spline configuration interaction method. Present values of the 
static dipole polarizabilities in the length, velocity and acceleration gauges are in good agreement 
with the best Hylleraas results. Also the tune-out wavelengths in the range from 400 nm to 4.2 gm 
for the four lowest triplet states are identified, and the magic wavelengths in the range from 460 nm 
to 3.5 gm for the 2^5 —)■ 3^S, 2^5 —>■ 2^P, and 2^5 —)■ 3^P transitions are determined. We show 
that the tune-out wavelength of 2^5' state is 413.038 28(3) nm, which corroborates the value of 
Mitroy and Tang (Phys. Rev. A 88, 052515 (2013)), and the magic wavelength around 1066 nm 
for the 2^5 —>■ 3^P transition can be expected for precision measurement to determine the ratio 
of transition matrix elements (2^5 —)■ 2^P)/{3^P —)■ 6^5). 
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I. INTRODUCTION 


Precise calculations of dynamic dipole polarizabilities for atoms are of interest due to its 
importance in a number of applications. First, dynamic dipole polarizabilities can be used 
directly to analyse the ac Stark shift to pursue higher-precision atomic clocks [l|, . Second, 

investigation of the dynamic dipole polarizabilities can derive the magic wavelengths and 
tune-out wavelengths, which open a new route to determine the line strength ratio @,i 
and to test the relativistic and quantum electrodynamic (QED) effects upon the transition 
matrix element not on the energy 0,0. And at last, since both of the trapping potential 
depth and the photon-scattering rate are dependent on the polarizabilities, the calculations 
of the dynamic dipole polarizabilities can provide reliable reference for experimental design 
to trap atoms in efficiency 0,0. 

As the simplest two-electron system, the accurate theoretical calculations and experi¬ 
mental measurement^ of the energy levels for helium can be used to test the the three-body 
bound QED theory 


12 


to determine the fine structure constant with high- 
to extract the nuclear information without resortin g to any model 
to develop the multi-electron atomic structure theory 0 _ 
transition 2^S 2^P and the doubly-forbidden transition 


precision 


10 - 


14 , and 


15 


Recently, the resonance 
2^S —)■ 2^5 of helium iso¬ 


topes have attracted great interest for the determination of nuclear charge radius differ- 
0,0 [l3. [l0 . Combined the laser cooling with magneto-optical trap techniques, the 


ence 


transitions 2^S —)■ 2and 2—)■ 3^P of helium are also demonstrated to produce high 
density quantum gas 17|] . The key point to improve the experimental measurement precision 


for helium is setting the laser frequency at the magic wavelength to eliminate effectively the 
ac Stark shift induced by the trap light. 

At present, there are lots of literatures focused on the accurate calculations of the energy 
and polarizabilities [l8-23| for the ground state of helium. For example, the non-relativistic 
ground-state energy has been achieved up to 46 digits [1^, and the static dipole polarizabil¬ 
ity of the ground-state helium, which includes the effect of mass polarization, the relativistic 
and leading QED corrections, has been determined to 1.383 191(2) within 2 ppm accu¬ 
racy [1^. However, compared with the ground state, there are very few calculations of 
dynamic polarizabilities for the triplet states of helium. As we known, for the metastable 
state 2 of helium, Glover et al. listed the rigorous upper and lower bounds of the dynamic 
dipole polarizabilities [l^. Chung provided dynamic polarizabilities for frequencies up to 


the second excitation threshold by using a variation-perturbation scheme [2^. Chen used 
a conhguration interaction (Cl) scheme with B-^line functions j29| to improve the conver¬ 


gence of the dynamic dipole polarizabilities |30,1^ . And Rerat et al. presented the dynamic 


dipole polarizabilities of helium at both real and imaginary frequencies using time-dependent 


gau ge-in variant method 
ties 


_|. In 2005, Chernov et al. calculated the dynamic polarizabili- 

33[ by using the quantum defect Green function formalism. For others triplet states of 


helium, there are fewer reports can be referenced 33|, |34 


In this work, hrstly, we have performed the calculations of static dipole polarizabilities 
for the low-lying triplet states 2^S, 2^P, 3^S and 3^P of helium with the conhguration 
interaction method based on B-spline functions in the length, velocity and acceleration 
gauges. Then the dynamic dipole polarizabilities of 2 ^S, 2 ^P for frequencies below the 
second excitation threshold, and 3 ^S, 3 ^P for frequencies below the hrst ionization threshold 
are calculated utilizing oscillator strengths and energy differences obtained in the length 
gauge. In addition, using the dynamic dipole polarizabilities, the magic wavelengths for the 
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three transitions 2^5' —)■ 3^S, 2^S ^ 2^P, and 2^5" —?• 3^P, and the tune-out wavelengths 
for the four lowest triplet states 2 ^S, 2 ^P, 3 ^S, and 3 ^P are determined with high accuracy. 


II. DIPOLE POLARIZABILITY 


The dynamic dipole polarizability of the magnetic sub-level \LgMg) is 

^ ^ ^ ^ - Lg{Lg -M) 

= aiiu) + - ^ -«! (w), 


( 1 ) 


where ai{u;) and aj (w) are the dynamic scalar and tensor dipole polarizabilities respectively, 
they are expressed as the summation of all allowed-transition intermediate states, including 
the continuum. 


f (1) 
Jgn 


/ \ Jg 




( 2 ) 




L.+L, /30(24+1)4(24-1) / 1 1 2 


(2-^9 + 3)(-^3 + 1) \ ^g ^g J — uj"^ 


/i 


( 1 ) 

gn 


(3) 


In the above formula, AEgn is the transition energy between the initial state \NgLgMg) 
and the intermediate state |iV„L„M„), uj is the photon energy of external electric field, 
and the dipole oscillator strength fgn have different expressions in the length, velocity, and 
acceleration gauges respectively. 


f(i) = 
J gn 


‘^\{NgLg\\ E nC^^\h)\\NnLn)?AE,^ 


i=l,2 


gn 


3{2Lg + 1) 


(4) 


d 


^{NgLgW Y. ^C^^\h) + h{^u]^t)rYC^^\mNnLn)\\AEg^)-^ 

/-i \ i—\ 2 dv i 

-. (5) 


f(l) ^ 

J gn 


3(24 + 1) 

2|(iVAII E Z4'Ct‘)(fj)||lV„4)P(A4„)-3 


i=l,2 


3(24 +1) ’ 

where Y being the electronic dipole transition operator of two-electron system, 

i=l,2 

Z is the nuclear charge number, or are the orbital quantum number of a electron, and 
the function b{i; A) is dehned as 


= e + 1 

b{i-i+l) = -1. 


(7) 


According to the Eqs. ([2]) and ([3]), for the case of = 0 initial state, the dynamic scalar 
and tensor dipole polarizabilities are 


ai{u) = ai{P,u) 
af(a;) = 0, 


( 8 ) 

(9) 
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where ai{P,uj) represents the contributions of the intermediate state with the angular mo¬ 
mentum number = 1. 

For the initial state of = 1, the dynamic scalar and tensor dipole polarizabilities are 
expressed as 


0 : 1 ( 0 ;) = oi(S', 0 ;)-|-oi(P, 0 ;)-|-oi(Zi), 0 ;) (10) 

al{u) = -ai{S,uj) + ]^ai{P,uj)-^ai{D,uj), ( 11 ) 

where oi(5', 0 ;) and ai(Zl,o;) are the contributions of the natural parity state {ss)S and 
{sd)D respectively, and oi(P, 0 ;) is the contribution of the unnatural parity state of {pp')P 
electron conhguration. 

In order to calculate the dynamic dipole polarizabilities, the fundamental atomic structure 
information of energies and wavefunctions are needed to obtain hrstly. In our calculations, 
the conhguration interaction method based on B-spline functions are adopted to get the 
energies and wavefunctions for helium. 


III. CONFIGURATION INTERACTION WITH B-SPLINE BASIS 


The Hamiltonian for two-electron system is given in second-quantized form as 

1 




— / Oji 




( 12 ) 


ijk£ 


where e* is the ith energy eigenvalue of the single-particle Schrddinger equation, gijki is 
two-particle matrix element of the Coulomb interaction, and and are creation and 
annihilation operators for the ith electron respectively. The single-particle quantum state is 
presented as \niiimimsi), here n* is the principal quantum number, ii is the orbital angular 
momentum, m* and rusi ^he orbital and spin angular momentum projection, respectively. 

The two-electron wavefunction 'ipij{LS) is expressed as a linear combination of 
conhguration-state wavefunctions (pij^LS), 


= , (13) 

ij 

and the conhguration-state wavefunction has the following expression. 


0ij(f^*S*) TJij 


E 


E 


{iiTiii, ijmj\LM){l/2msi-, l/2ms^\SMs)alapO ), (14) 


rriimj ms-iris- 


where rjij is a normalization constant given by 


V^j 


1, * ^ J 

= i. 


(15) 


The Clebsch-Gordan coefficients {iimf, ijmj\LM) and {l/2msi',l/2msj\SMs) represent ii 
and ss coupling, respectively, |0) is the vacuum state and aj^|0) represents the ith eigen- 
wavefunction of the single-particle Schrddinger equation with energy eigenvalue e^. The 
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configuration-state wavefunctions are independent of magnetic quantum numbers of mj, mj, 
rUsi - From the interchange symmetery of the Clebsch-Gordan coefficients, it follows 

that 

, (16) 

which implies (j)ii{LM) = 0 unless L + S is even. 

According to the expansion of the wavefunctions, the matrix elements of Hamiltonian is 

(f/’ij (FAf) Ihf ^ A ^ ^ ^ 

ki ij,k£ 

where the potential energy matrix element Vij^ke between different configurations is 


^ijki ^ij^kl 


;; ei t 




The quantity Xy{ijktj in the above equation is given by 

where (£i||C"^||£fc) is angular reduced matrix element, 

{(iWC-Wk) = (-1)'‘V(24 + 1)(24 + 1) I 0 0 0 


(18) 


(19) 


( 20 ) 


The two-electron radial integral R,^{ijki) of the Coulomb interaction is written as 

R^{ijki) = J J rfdrirldr2Ri{ri)Rk{ri)-^Rj{r2)Reir2) , ( 21 ) 

where r< and r> are the minimum and maximum of ri and r 2 , and Ri{r) is the radial 
wavefunction of the fth single-electron orbital. 

Using the variational method, the followed configuration interaction equations can be 
obtained, 

^ ^ T ^j)^ik^jl T ^ijki\(^kl 
kt 

where A and Cij are the eigen-energy and eigen-wavefunction for two-electron atoms, respec¬ 
tively. 

Before solved the Cl equations, the energies and wavefunctions for single-electron orbital 
are obtained firstly. In our calculations, B-splines are used to expand the radial wavefunction 
for the ith single-electron orbital. 


R,(r) = 5^c;B,(r). 

j 


(23) 
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TABLE I: Convergence of the energies (in a.u.) for the metastable state 2^S and the oscillator 
strengths ™ length gauge of helium as the number of B-splines N and partial waves 

(■max increased. The numbers in parentheses of the extrapolated values give the computational 


uncertainties. 




Energy 



f(L) 

■>2 3 S^2^P 


^max 

N=30 

N=35 

N=40 

N=30 

N=35 

N=40 

2 

-2.175 220 4147 

-2.175 220 4306 

-2.175 220 4345 

0.539 818 2380 

0.539 818 2108 

0.539 818 2056 

3 

-2.175 227 0950 

-2.175 227 1164 

-2.175 227 1220 

0.539 204 6079 

0.539 204 5519 

0.539 204 5388 

4 

-2.175 228 5828 

-2.175 228 6093 

-2.175 228 6165 

0.539 117 1546 

0.539 117 0681 

0.539 117 0460 

5 

-2.175 229 0255 

-2.175 229 0559 

-2.175 229 0647 

0.539 096 9772 

0.539 096 8626 

0.539 096 8313 

6 

-2.175 229 1847 

-2.175 229 2183 

-2.175 229 2283 

0.539 090 7709 

0.539 090 6327 

0.539 090 5931 

7 

-2.175 229 2501 

-2.175 229 2859 

-2.175 229 2970 

0.539 088 4620 

0.539 088 3052 

0.539 088 2583 

8 

-2.175 229 2796 

-2.175 229 3171 

-2.175 229 3289 

0.539 087 4793 

0.539 087 3088 

0.539 087 2559 

9 

-2.175 229 2939 

-2.175 229 3326 

-2.175 229 3450 

0.539 087 0173 

0.539 086 8365 

0.539 086 7790 

10 

-2.175 229 3013 

-2.175 229 3407 

-2.175 229 3536 

0.539 086 7826 

0.539 086 5947 

0.539 086 5336 

Extrap. 


-2.175 229 36(2) 



0.539 086 4(3) 



where {c*} are the expansion coefficients, and the following exponential knots are employed, 


' ti = 0, 

(i+k—l 


exp 
Rq - 


(tRi) 

exp[7i?o] - 1 


- 1 


5 


i = 

i = N + k — l,N + k. 


( 24 ) 


where Rq is the box size, which need to be chosen large enongh to make sure the contributions 
to dynamic dipole polarizabilities from higher excited-state are included, especially when the 
photon energy cj is large. The non-linear parameter 7 is also need to be adjusted to get 
more accurate ground-state energy of helium, then the value of 7 is fixed the same for all 
the triplet states to simplify the integral of B-splines. 


IV. RESULTS AND DISCUSSIONS 

In our calculation, Rq = 200 a.u. and j = Rq x 0.038 are used throughout the paper. 
Using the fixed values of Rq and 7 , we get the ground-state energy —2.879 0284 a.u. under 
S-wave approximation with 30 B-splines of order 7, which has 7 significant digits with the 
S-wave limit value —2.879 028 767 29 a.u. jUj. 


A. Energies and Oscillator Strengths 

Table [I] is the convergence of the energies for the metastable state 2^S and the oscillator 
strengths of 2 —)■ 2 transition in the length gauge of helium as the number of basis set 

and partial waves increased. For the energy, increase of the number of partial wave change 
less than the number of B-spline N increased. This convergent style for the energy suggested 
that we can fix partial wave (in our work we fix (max = 10 ), then increase the number of 
B-spline N to avoid too enormous number of CL Considering both the effect from N and 
(max, the extrapolated values are given in the last line of the Table [H The hnal converged 
value for the energy is —2.175 229 36(2) a.u., which is in excellent agreement with the result 
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TABLE II: Comparison of the energies (in a.u.) in the lengthy gauge for the four lowest triplet 
states of helium. The numbers in the parentheses are the computational uncertainties. 


State 

Present 

Ref. [3^ 

Ref. [3^ 

Ref. 

2^S 

-2.175 229 36(2) 

-2.175 228 8 

-2.175 229 378 176 

-2.175 229 378 236 791 30 

3^S 

-2.068 689 07(2) 

-2.068 688 8 

-2.068 689 067 283 

-2.068 689 067 472 457 19 

2^P 

-2.133 164 17(2) 

-2.133 163 4 

-2.133 164 181 6 

-2.133 164 190 779 273(5) 

3^P 

-2.058 081 08(2) 

-2.058 080 6 

-2.058 081 077 2 

-2.058 081 084 274 28(4) 


TABLE III: Comparison of the oscillator strengths in three different gauges for helium. The 


numbers in the parentheses are the computational uncertainties. The values of Ref. 3^ is the 


average of the length and velocity gauges, the results of Ref. [37| are in the length gauge, and 


Ref. 


39l| gives the values in three different gauges. 




Present 


Ref. [361 

Ref. [37] 


Ref. [3^ 


Transition 

f[^) 

Jgn 

Tun 

Jgn 

Jgn 



Jgn 

A'^ ) 

Jgn 

Jgn 

2^S 

0.539 0864(3) 

0.539 0865(2) 

0.539 078(6) 

0.5391 

0.539 0861 

0.5392(8) 

0.539(3) 

0.56(3) 


0.890 8518(2) 

0.890 8518(4) 

0.890 83(3) 

0.8910 

0.890 8513 

0.890(2) 

0.889(7) 

0.85(6) 

2^P ^ 3^D 

23p^ 2 3pe 

0.610 2255(2) 
0.180 480 28(2) 

0.610 2255(2) 
0.180 4803(2) 

0.610 2247(3) 
0.180 4803(3) 

0.610 24 

0.610 2252 

0.611(2) 

0.609(2) 

0.609(3) 

33p^ 4^D 
3^P^ 

0.477 5943(2) 
0.135 420 99(2) 

0.477 5943(2) 
0.135 420 99(3) 

0.477 593(3) 
0.135 420 98(4) 

0.477 60 

0.477 5938 

0.474(3) 

0.476(1) 

0.494(5) 


-2.175 229 378 176 a.u. of Cann and Thakkar [36|. The extrapolated oscillator strength 


0.539 0864(3) has 6 signihcant digits with the value 0.539 0861 of Drake 

A similar convergence pattern exists for the energies and oscillator strengths in the veloc¬ 
ity and acceleration gauges for the other triplet states of helium. The hnal convergent results 
of the energies in the len gthy gauge are presented in Table [TTl Our energies are much more 


jtnv 

accurate than the values 38| by two order of magnitudes, which are obtained by using the 
B-splines Cl method with different number of conhguration states. And our results for the 
2 and 2 states have 8 signihcant digits with the explicitly correlated calculations 36 
and the Hylleraas results [^ . 

Table IIIII lists the comparison of the oscillator strengths for some selected transitions. 
The superscripts (L), (R), and (A) represent results obtained in the length, velocity and 
acceleration gauges, respectively. For the dipole oscillator strength of 2 —)■ 2 transition, 
the value in the acceleration gauge is less accurate than the results from length and velocity 
gauges, but our results for 2^S —)■ 2^P in three gauges are correspondingly much more 


accurate than the values in different gauges of Ref. [39| by three order of magnitudes. All of 


our results in Table m are much more accurate than the previous values of Refs. |36|, 


And for the oscillator strengths of other transitions, our results in the length and velocity 
gauge are in excellent agreement with the Hylleraas calculations of Drake In addition, 
the oscillator strengths from the initial states 2 ^P and 3 ^P transit to the unnatural parity 
states 2 ^P® and 3 ^P® are also listed in the Table IIIII 
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TABLE IV: Convergence of the static dipole polarizabilities (in a.u.) in three different gauge for 
the metastable state 2 of helium as the number of B-splines N and partial waves imax increased. 


The numbers in parentheses give the computational uncertainties. 





aW(0) 










(0) 



^max 

N=30 

N=35 

N=40 

N=30 

N=35 

N=40 

N=30 

N=35 

N=40 

2 

315.433 

490 

315.433 

397 

315.433 

373 

315.171 

086 

315.170 

980 

315.170 

949 

312.742 

586 

312.745 

210 

312.745 

924 

3 

315.606 

026 

315.605 

905 

315.605 

872 

315.571 

449 

315.571 

315 

315.571 

274 

315.072 

959 

315.076 

939 

315.078 

064 

4 

315.626 

281 

315.626 

136 

315.626 

095 

315.618 

705 

315.618 

550 

315.618 

502 

315.463 

517 

315.468 

762 

315.470 

301 

5 

315.630 

187 

315.630 

024 

315.629 

976 

315.627 

938 

315.627 

769 

315.627 

715 

315.562 

445 

315.568 

761 

315.570 

684 

6 

315.631 

213 

315.631 

037 

315.630 

984 

315.630 

396 

315.630 

219 

315.630 

160 

315.594 

442 

315.601 

601 

315.603 

855 

7 

315.631 

548 

315.631 

364 

315.631 

306 

315.631 

205 

315.631 

023 

315.630 

960 

315.606 

534 

315.614 

322 

315.616 

846 

8 

315.631 

678 

315.631 

487 

315.631 

426 

315.631 

515 

315.631 

330 

315.631 

265 

315.611 

603 

315.619 

840 

315.622 

577 

9 

315.631 

735 

315.631 

540 

315.631 

476 

315.631 

649 

315.631 

462 

315.631 

396 

315.613 

883 

315.622 

434 

315.625 

333 

10 

315.631 

763 

315.631 

565 

315.631 

500 

315.631 

712 

315.631 

524 

315.631 

457 

315.614 

962 

315.623 

727 

315.626 

746 

Extrap. 



315.631 

5(2) 





315.631 

4(2) 





315.63(2) 




TABLE V: Comparison of the static dipole polarizabilities (in a.u.) for helium. The numbers in 
parentheses give the computational uncertainties. 


state 

a[^\0) 

aS'bo) 


Ref. [40] 

2^S 

315.6315(2) 

315.6314(2) 

315.63(2) 

315.631 47(1) 

3^S 

7937.584(2) 

7937.583(2) 

7937.4(2) 

7937.58(1) 

2^P 

46.70793(4) 

46.70794(4) 

46.71(2) 

46.707 7482(3) 

3^P 

17305.67(3) 

17305.67(4) 

17311(2) 

17305.598(3) 


B. Static Dipole Polarizabilities 


Table IIVI gives the convergence of the static dipole polarizabilities for the metastable 
state 2 of helium as the number of basis set and partial waves increased, and the last line 
lists the extrapolated values. From this table, we can see in the length and velocity gauges, 
the convergence style are the same, the results are decreased as the number of basis sets 
N increased for a same imax- However in the acceleration gauge, the values are increased 
as the number of basis sets N increased for a same imax- The hnal convergent value in the 
length gauge is 315.631 5(2), which has 6 signihcant digits compared with the most accurate 


Hylleraas value 315.631 47(1) of Yan |40 


Table |V] gives the comparison of static dipole polarizabilities for the four lowest triplet 
states of helium. The results between the length and velocity gauges are in perfect agreement. 
The values obtained in the acceleration gauge are less accurate than the results of lengthy 
and velocity gauges by two order of magnitudes. Present results for the 2 and 3 states 
in the length and velocity gauges agree with the Hylleraas values at the 10“^ level, and 
our values for the 2 and 3 states in the length and velocity gauges agree with the 
Hylleraas values at the 10“® level. For the acceleration gauge, present for 2^P 

and 3^P states just have 3 singihcant digits compared with the Hylleraas values |40|. 


C. Dynamic Dipole Polarizabilities 

Table |VT] lists the dynamic dipole polarizabilitity for the metastable state 2 of helium 
for some selective frequency from 0 to 0.12 a.u., the hgures in parentheses represent com- 
















putational uncertainties. It seen clearly from this table, all of our values have at least 5 
significant digits except the results of a; = 0.04 a.u., uj = 0.110 a.u., and uj = 0.115 a.u., 
which only have 4 significant digits. That’s because there is always a tune-out wavelength lo¬ 
cated in the vicinity of these positions j^, the relativistic and finite nuclear mass corrections 
may effect the uncertainities of dynamic dipole polarizabilities. 

Table I VII al so makes a comparison of the present results with available values from other 


literatures [28|, [Sll l32|]. All of ours results lie within the boundary of Glover and Wein- 
hold’s ji^, which gives the rigorous upper and lower limits for the dynamic dipole polar¬ 
izability at a wide frequency range. In the low-frequency region, our values are in good 


agreement with Ref. [SlJ, which are also obtained by using B-spline Cl method. For exam¬ 


ple, ours values have the same hve signihcant digits as theirs. As the frequency u increased, 
the differences between present results and values of Ref. jsij increased, especially for the 
u = 0.115 a.u., the difference of the dynamic polarizabilities can reach to about 2.2 Og. The 
reason for this is that the box size Rq of B-spline adopted in present and Ref. j3l| calculations 
are different. As u increased, Rq should be chosen big enough to make sure the transition 
to high-excited states, especially the transition energies of those excited states near ca, can 
be included in the calculation of polarizabilities. For example, if the box size Rq = 200 a.u. 
is adopted, we get ai(0.140) = 1.91012(92) Og and ai(0.145) = —65.327(5) Og. If the box 
size is set as Rq = 50 a.u., then we get q;i( 0.140) = 1.57(6) aj] and ai(0.145) = —68.6(6) Og, 
which are less accurate than the values from Rq = 200 a.u.. However, the box size is not 
the bigger the better for the B-spline Cl calculation, oppositely, the loss of accuracy will 
occur under the same number of B-spline for the big box size. In order to get more accurate 
values, the number of B-spline should be increased, which makes the number of Cl increased 
exponentially and slows down the convergent process of our calculations. So in our practical 
calculation, we need to chose appropriate Rq to get accurate value for large oj and to avoid 
large number of Cl at the same time. 

Table IVIII lists some selective values of dynamic dipole polarizabilities of 2 ^P, 3 ^S, 
and 3 states for the He atom. For the 2 ^P state, we calculate the dynamic dipole 
polarizabilities for frequency below the second excitation threshold, and for the 3^5 and 
3 ^P states, we only list the dynamic dipole polarizabilities for frequency u below the hrst 
ionization threshold. All of our results are very accurate except few values for the frequency 
near resonance transition energy or ionization threshold. 

The dynamic dipole polarizabilities for the lowest four triplet states of helium are also 
plotted in the Figs. [1]-[6] as the photon energy u. For the non-zero angular momentum state, 
the polarizability depends upon its magnetic quantum number M because of both scalar 
and tensor polarizabilities existing, so the dynamic dipole polarizabilities for 2 and 3 ^P 
states are divided into two cases as M = 0 and \M\ = 1. The crossing points between 
a curve and the horizontal zero line are called as tune-out wavelengths, which denoted as 
solid magenta circle, and the crossing points between two curves are the magic wavelengths, 
which denoted as blank red circle. The vertical lines are the resonance transition positions. 


D. Tune-Out Wavelengths 

Table IVIHI lists the values of tune-out wavelengths in the 400-4200 nm region, which 
marked as solid magenta circle in the Figs. [U-El For the metastable state of helium, Mitroy 
and Tang have obtained the 413.02(9) nm tune-out wavelength by incorporating Hyller- 
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TABLE VI: Comparison of the dynamic dipole polarizabilities (in a.u.) of the 2^S state for the 
He atom. The numbers in parentheses give the computational uncertainties. 


Ijj 

Present 

Ref. [27] 

Ref. [28] 

Ref.[^ 

Ref. [3^ 

0.000 

315.6315(2) 

(315.61, 316.83) 

315.63 

315.630 

315.92 

0.005 

320.0105(2) 

(319.99, 321.23) 

320.01 

320.009 

320.31 

0.010 

333.9323(2) 

(333.91, 335.21) 

333.93 

333.931 

334.25 

0.015 

360.1322(2) 

(360.10, 361.53) 

360.12 

360.130 

360.50 

0.020 

404.8286(2) 

(404.79, 406.43) 

404.81 

404.825 

405.28 

0.025 

482.3411(3) 

(482.29, 484.29) 

482.31 

482.335 

482.95 

0.030 

631.4758(3) 

(631.39, 634.10) 

631.42 

631.463 

632.41 

0.035 

1001.751(2) 

(1001.53, 1006.08) 

1001.59 

1001.71 

1003.68 

0.040 

3192.7(2) 

(3190.16, 3207.18) 

3190.67 

3192.17 

3205.28 

0.045 

-2097.602(3) 

(-14717.51, -2050.70) 

-2098.66 

-2097.89 

-2097.0 

0.050 

-725.4477(2) 

(-729.45, -718.73) 

-725.60 

-727.490 

-726.27 

0.055 

-416.4749(2) 

(-419.07, -413.67) 

-416.54 

-416.492 

-417.20 

0.060 

-281.12719(2) 

(-283.14, -279.41) 

-281.16 

-281.137 

-281.78 

0.065 

-205.59918(2) 

(-207.31, -204.30) 

-205.63 

-205.606 

-206.21 

0.070 

-157.5702(2) 

(-159.14, -156.44) 

-157.59 

-157.575 

-158.17 

0.075 

-124.33628(2) 

(-125.87, -123.25) 

-124.35 

-124.340 

-124.94 

0.080 

-99.86344(2) 

(-101.47, -98.72) 

-99.88 

-99.867 

-100.49 

0.085 

-80.8752(2) 

(-82.67, -79.60) 

-80.89 

-80.878 

-81.55 

0.090 

-65.36352(2) 

(-67.54, -63.84) 

-65.37 

-65.366 

-66.11 

0.095 

-51.88420(2) 

(-54.81, -49.93) 

-51.92 

-51.888 

-52.74 

0.100 

-39.05227(2) 

(-43.49, -36.31) 

-39.11 

-39.059 

-40.07 

0.105 

-24.69053(2) 

(-32.65, -20.33) 

-24.82 

-24.709 

-26.03 

0.110 

-2.1515(2) 


-2.66 

-2.248 


0.115 

93.381(2) 

(-6.27, 128.35) 

84.66 

91.175 

74.32 

0.120 

-125.6637(4) 


-137.94 




aas matrix elements for the transition to 2^P and 3^P manifolds and core-polarization 
model matrix elements for other transitions, and they predicted the tnne-ont wavelength 
aronnd 413 nm can be nsed to test the QED effect. Recently, a experimental measnrement 
of Ken Baldwin’s gronp report the tnne-ont wavelength being 413.0938(9Stat. )(20Syst.) 
nm and another theoretical calculation by Notermans et al gives 414.197 nm tune-out 
wavelength by using available tables of level energies and Einstein A coefficients. Our tnne- 
ont wavelength of ab-initio calculation is 413.038 28(3) nm, which corroborates the value 
413.02(9) nm of Mitroy and Tang [^. The difference between the theoretical calculations 
and the experimental measurement may caused by hnite nuclear mass, relativistic and QED 
corrections, which calls for great efforts for theoretical calculation to improve the precision 
for QED test. 
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TABLE VII: Dynamic polarizabilities (in a.u.) of 3^S, 2^P, and states for the He atom. The 
numbers in parentheses give the computational uncertainties. 


state 

335 

2 

3p 


3^P 

ijj 

ai(w) 

ai(w) 

(w) 

ai{u)) 

ai\uj) 

0.000 

7937.584(2) 

46.70793(3) 

69.5964(2) 

17305.67(3) 

336.768(3) 

0.005 

10199.363(2) 

45.83274(4) 

70.8984(2) 

-8051.87(2) 

3534.464(2) 

0.010 

71124.62(5) 

42.97603(4) 

75.0559(2) 

-23498.47(2) 

23293.67(2) 

0.015 

-7892.735(2) 

37.3217(2) 

82.9492(2) 

3942.877(2) 

-3345.9627(4) 

0.020 

-3062.782(2) 

26.972(2) 

96.59195(3) 

4749.082(4) 

-3445.013(2) 

0.025 

-1685.264(2) 

7.4367(2) 

120.65874(3) 

4962.08(2) 

-246.35(2) 

0.030 

-1031.851(2) 

-33.77638(2) 

167.92335(6) 

-2249.18(2) 

111.731(2) 

0.035 

-305.1046(2) 

-145.8045(2) 

287.9533(2) 

502.52(2) 

-765.454(5) 

0.040 

-792.8559(2) 

-860.666(3) 

1013.44(2) 

-1078.384(2) 

53.0407(4) 

0.045 

-512.4321(2) 

924.3286(4) 

-757.2164(4) 

-1078.7(2) 

78.96(2) 

0.050 

-602.5701(5) 

498.7080(2) 

-311.6507(2) 

152(2) 

14.1(2) 

0.055 

-703.68(2) 

447.8331(3) 

-231.6643(2) 



0.060 

-620.32(3) 

517.4623(6) 

-255.4211(4) 



0.065 


-563.39(2) 

907.7(2) 



0.070 


561.824(2) 

-27.2995(3) 



0.075 


1635.892(8) 

-167.107(2) 



0.080 


-1497.409(7) 

137.702(2) 



0.085 


-418.409(2) 

26.1352(2) 



0.090 


-185.5369(3) 

-1.0091(2) 



0.095 


-29.424(2) 

-39.1638(3) 



0.100 


211.866(3) 

-18.9742(3) 





0.035 0.040 0.045 0.050 0.055 0.060 


photon energy (a.u.) 

FIG. 1: (Color online) Dynamic dipole polarizabilities (in a.u.) of helium for the photon energy 
0.035 < CO < 0.06 a.u. The solid black line denotes the dynamic polarizabilities for 2^S state, and 
the dashed blue line represents the dynamic polarizabilities for 3^5 state. The crossing points 
denoted as solid magenta circle are the tune-out wavelengths, and the crossing points marked as 
blank red circle are the magic wavelengths. The vertical lines are the resonance transition positions, 
and the green line is a horizontal zero line. 
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FIG. 2: (Color online) Dynamic dipole polarizabilities (in a.u.) of helium for the photon energy 
0.06 < w < 0.10 a.u. The solid black line denotes the dynamic polarizabilities for 2 state, and the 
dashed blue line represents the dynamic polarizabilities for 2^P{M = 0) state. The crossing points 
denoted as solid magenta circle are the tune-out wavelengths, and the crossing points marked as 
blank red circle are the magic wavelengths. The vertical lines are the resonance transition positions, 
and the green line is a horizontal zero line. 



photon energy (a.u.) 


FIG. 3: (Color online) Dynamic dipole polarizabilities (in a.u.) of helium for the photon energy 
0.06 < tu < 0.10 a.u. The solid black line denotes the dynamic polarizabilities for 2^S state, 
and the dashed blue line represents the dynamic polarizabilities for 2^P(|M| = 1) state. The 
crossing points denoted as solid magenta circle are the tune-out wavelengths, and the crossing 
points marked as blank red circle are the magic wavelengths. The vertical lines are the resonance 
transition positions, and the green line is a horizontal zero line. 


E. Magic Wavelengths 


The magic wavelength is the wavelength at which the polarizability difference for a transi¬ 
tion goes to zero, which means the hrst-order Stark shifts for the npper and lower levels of a 

Table HXl presents all the valnes of magic wavelengths in the 


transition are the same 42 


460-3500 nm region marked in blank red circle in the Figs. [T]l6l The corresponding dynamic 
dipole polarizabilities at the magic wavelengths are also given in the last colnmn. For the 
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FIG. 4: (Color online) Dynamic dipole polarizabilities (in a.u.) of helium for the photon energy 
0.02 < w < 0.04 a.u. The solid black line denotes the dynamic polarizabilities for 2 state, and the 
dashed blue line represents the dynamic polarizabilities for 3^P(M = 0) state. The crossing points 
denoted as solid magenta circle are the tune-out wavelengths, and the crossing points marked as 
blank red circle are the magic wavelengths. The vertical lines are the resonance transition positions, 
and the green line is a horizontal zero line. 



FIG. 5: (Color online) Dynamic dipole polarizabilities (in a.u.) of helium for the photon energy 
0.04 < uj < 0.05 a.u. The solid black line denotes the dynamic polarizabilities for 2 state, and the 
dashed blue line represents the dynamic polarizabilities for 3^P(M = 0) state. The crossing points 
denoted as solid magenta circle are the tune-out wavelengths, and the crossing points marked as 
blank red circle are the magic wavelengths. The vertical lines are the resonance transition positions, 
and the green line is a horizontal zero line. 


magic wavelength of 1066.197(2) nm, there exist two terms, which play major contribution 
of the dynamic dipole polarizabilities for the 2^S and 3^P{M = 0) states respectively. Ta¬ 
ble |X]lists some contributions from different intermediate states for the 2 —>■ 3 ^P{M = 0) 
transition in detail at the magic wavelength of 1066.197(2) nm. We can see that the con¬ 
tribution from 2^P state to the polairzability of 2^S is about 99.87%, and the contribution 
from 6^5* state to the polairzability of 3^P{M = 0) is about 98.87%. According to the 
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FIG. 6: (Color online) Dynamic dipole polarizabilities (in a.u.) of helium for the photon energy 
0.01 < CO < 0.05 a.u. The solid black line denotes the dynamic polarizabilities for 2^S state, 
and the dashed blue line represents the dynamic polarizabilities for 3^P(|M| = 1) state. The 
crossing points denoted as solid magenta circle are the tune-out wavelengths, and the crossing 
points marked as blank red circle are the magic wavelengths. The vertical lines are the resonance 
transition positions, and the green line is a horizontal zero line. 


definition of magic wavelength a23s{oJm) = Oi^sp^com), we have the expanded form. 


/2 s^2 p + a2 3 s{Remainder;Um) = ^ —— + a^sp{Remainder 


^^23s^23p 




(25) 


where the second term in the left of Eq. fl25ll is all the contributions from other n^P states 
to the dynamic dipole polarizability of 2^S state, and the second term in the right of 
Eq. fl25p is all the contributions from other n^S, n^D, and n^P^ states to the dynamic 
dipole polarizability of 3^P{M = 0) state. If all the remainder terms are neglected, then 
the ratios of oscillator strengths and reduced matrix elements are written as 


f23S^23p _ 3(AE^3g^^3p uj^) 

fs^P-^-G^S ^^33p^(}3s 

|(23P|| ^ r,GW(f,)||23p)| 
M 2 35^2 3p _ 1=1,2 

~ |(33P|| ^ r,G(b(f,)||63P)| 

i=l,2 


(26) 


AP 3 3p^g3g( AP 2 35'^2 3p ^m) 

Ap2 35^2 3p(APg3p^g3^ 


(27) 


Combined present energy difference and the magic wavelength 1066.197(2) nm, the ratios 
of the oscillator strengths and the reduced matrix elements are determined and listed in 
Table IXII Present^ are the values of our ab-initio calculation, and Present^ are derived 
by substituting our theoretical energies and the magic wavelength of 1066.197(2) nm into 
the Eqs. fEB]) and (127)) . Compared with the explicitly correlated results of Ref. [s^, we 
believe our values of Present^ are reliable, since present oscillator strengths for 2^5* —)■ 2^P 
and 3^P —)■ transitions are much more accurate than the values of Ref. 36(] by at 
least one order of magnitude. In order to test the accuracy of the values derived from 
Eqs. fl2^ and fl27)) . we can compare the results between Present^ and Present^. It’s clearly 
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TABLE VIII: Tune-out wavelengths (At) for 2^S, 3^S, 2^P, and 3^P states of helium. The second 
and third column are the tune-out wavelengths in atomic units and nanometer respectively. The 
numbers in parentheses give the computational uncertainties. 


State 

a;t(a.u.) 

At(nm) 


0.110 312 66 ( 2 ) 

413.038 28 ( 3 ) 

3^,5 

0.035 488 102(2) 

1283.905 03(4) 


0.047 681 245(3) 

955.582 28(7) 


0.054 260 56(3) 

839.714(2) 


0.058 1756(2) 

783.204(2) 

2^P{M = 0) 

0.066 652 71(2) 

683.5934(2) 


0.094 341 03(2) 

482.9644(2) 


0.098 338 54(2) 

463.3316(2) 

2^P{M = ±1) 

0.097 382 82(2) 

467.8788(2) 

3^P{M = 0) 

0.023 315 997(3) 

1954.1670(3) 


0.033 637 088(2) 

1354.5570(2) 


0.036 548 517(2) 

1246.6539(3) 


0.042 005 02(2) 

1084.712(2) 


0.043 390 63(3) 

1050.074(2) 


0.046 5939(2) 

977.883(2) 


0.047 4031(2) 

961.189(3) 


0.049 4427(2) 

921.539(4) 


0.049 9644(2) 

911.917(4) 

3^P{M = ±1) 

0.010 911 33(2) 

4175.783(4) 


0.035 787 67(2) 

1273.1580(2) 


0.043 149 27(2) 

1055.9472(3) 


0.047 2945(1) 

963.397(2) 


0.049 9054(1) 

912.995(3) 


seen that the derived values 64.6653 and 4.677847 from the Eqs. fl26l) and fl2711 are in good 
agreement with our ab-initio values 65.48(2) and 4.7073(3) at the level of 1.3% and 0.7% 
accuracy respectively. If increasing the number of B-spline basis sets, and also considered 
the contribution of the remainder term, then improvement of the accuracy for the transition 
matrix elements ratio {M23s^2^p )/up to 0.5% is achievable. 

As we known that, present experimental technique is very difficult to measure matrix 
elements accurately, only 1% accuracy for one or two of the lowest transitions have been 


reported [^, 1^. Recently, Herold et al. present a method for accurate determination of 


5s — 6p matrix elements in rubidium by measurements of the ac Stark shift around tune-out 
wavelength 0. In our calculation, the particular magic wavelength around 1066 nm can 
be used for experiment measurement to determine the atomic transition matrix elements 
involved highly excited states for helium. 
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TABLE IX: Magic wavelengths (Am) for 2^S 3^S, 2^S 2^P{M = 0,±1) and 2^S ^ 


3^P(M = 0,±1) transitions of the He atom. The corresponding dynamic dipole polarizabilities 
at the magic wavelengths are given in the last column. The numbers in parentheses give the 
computational uncertainties. 

Transition 

iJmia.U.) 

\m(nm) 

ai{ujm)(a.u.) 

2^S 3^S 

0.036 004 592(2) 

1265.48724(4) 

1151.058(2) 


0.048 775 162(5) 

934.1507(2) 

-872.007(2) 


0.058 0244(2) 

785.245(3) 

-324.31(2) 

2'^S 2^P(M = 0) 

0.066 228 54(2) 

687.9716(2) 

-191.9157(2) 


0.093 2295(2) 

488.7225(2) 

-56.5068(2) 


0.098 0166(2) 

464.8532(3) 

-44.169 53(4) 

2^S 2^P(M = ±1) 

0.095 7473(2) 

475.870 85(6) 

-49.962 90(4) 

2^S 3^P(M = 0) 

0.023 434 415(3) 

1944.2923(3) 

453.0297(2) 


0.034 618 094(3) 

1316.1716(2) 

955.639(2) 


0.037 202 6377(2) 

1224.73446(2) 

1410.365(2) 


0.042 734 44 ( 3 ) 

1066 . 197 ( 2 ) 

-9487(2) 


0.044 532 37(2) 

1023.151(2) 

-2511.07(3) 


0.047 069 3(2) 

968.007(3) 

-1196.32(3) 


0.048 478 5(2) 

939.867(4) 

-915.97(3) 


0.049 762 5(3) 

915.616(5) 

-750.19(3) 

2^S 3^P(M = ±1) 

0.013 031 86(2) 

3496.304(3) 

348.0648(2) 


0.037 297 86(2) 

1221.607 75(2) 

1436.592(2) 


0.044 421 15(2) 

1025.713(2) 

-2633.59(3) 


0.048 309 5(2) 

943.156(4) 

-942.89(2) 


V. CONCLUSIONS 

The calculations of the energies and the main oscillator strengths for the four triplet 
states (2 ^S, 3 ^S, 2 ^P, and 3 ^P) in the length, velocity and acceleration gauges are carried 
out by the configuration interaction based on the B-spline functions. Also the accurate 
dynamic dipole polarizabilities for the four lowest triplet states are obtained. Ours static 
dipole polarizabilities in the length and velocity gauges have 5-6 significant digits, which 
are in excellent agreement with the variational Hylleraas calculations. Present work lays 
solid foundation for the further to calculate the relativistic and QED effects on the dynamic 
polarizabilities of helium. 

In particular, the tune-out wavelengths for the four triplet states and magic wavelengths 
for the three transitions of 2 —>■ 3 ^S, 2^S ^ 2 ^P, and 2^S —?• 3^P are determined 
with high precision. Our tune-out wavelength 413.038 28(3) nm of the metastable state 
validate the value of Mitroy and Tang . And the magic wavelength around 1066 nm for 
2 ^5 —)■ 3 transition is proposed for experimental measurement to determine the ratio of 
the transition matrix elements (2^5 —)■ 2^P)/(3^P —)■ Q^S), this is a unique way to obtain 
accurate transition matrix element involved highly excited states. Also we expected that 
other tune-out wavelengths and magic wavelengths can provide theoretical reference for the 
precision-measurement experiment design in the future. 
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TABLE X: Contributions from some of intermediate states to the dynamic dipole polarizability of 
2^S and 3^P{M = 0) states at the magic wavelengths 1066.197(2) nm. 


uj{a.u.) 

0.042 734 44(3) 

X{nm) 

1066.197(2) 

Intermediate states 

2^S 

2^P 

- 9499.066 

3^P 

5.418 


1.386 

Others 

5.234 

Total 

- 9487.028 

Intermediate states 

3^P{M = 0) 

3^S 

519.840 

435 

-320.109 

5^S 

-118.888 


- 9379.600 

3^D 

-73.9017 


-517.059 


-395.758 


506.415 

Others 

292.0327 

Total 

- 9487.028 


TABLE XL Comparison of the ratios for the oscillator strengths {f23S-^23p)/{f3^P-^6^s) and the 
reduced matrix elements {M 23 p)/{^ 3 ^p^q^ s)■ Present^ are the value of our ab-initio calcu¬ 
lation, and Present^ are derived by substituting our theoretical energies and the magic wavelength 
of 1066.197(2) nm into the Eas.(l26p and (j27p. The numbers in parentheses give the computational 
uncertainties. 


(/235^23p)/(/33p^63s) 

(M 2 35^2 3p)/ 

Present^ 

65.48(2) 

4.7073(3) 

Present^ 

64.6653 

4.677847 

Ref. [3^ 

65.5308 

4.789739 
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